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With a combined analysis of data on Y(5S) —y hb(lP,2P)n + TY~ and T(5 S) —¥ B^B^n in an 
effective field theory approach, we determine resonance parameters of Zb states in two scenarios. In 
one scenario we assume that Z\, states are pure molecular states, while in the other one we assume 
that Zb states contain compact components. We find that the present data favor that there should 
be some compact components inside Z\] ] associated with the molecular components. By fitting the 
invariant mass spectra of Y(5S) — > hb(lP , 2P)n + n~ and Y(5 S) —> B^B*ir, we determine that the 
probability of finding the compact components in Zb states may be as large as about 40%. 

PACS numbers: 


I. INTRODUCTION 

Two charged bottomonium-like states 2^(10610) and 2^(10650)—denoted as Z b and Z ' b —were discovered 
by the Belle Collaboration in decays T(55) —> Y(nS)TT + n~ for n = 1,2, or 3 and T(5S) —i hb(mP)ir + ir~ 
for to = 1 or 2 HI- The masses and decay widths averaged over the five channels are mzu = 10607.2 ± 2.0 
MeV, F Zb = 18.4 ± 2.4 MeV, and m z > b = 10652.2 ± 1.5 MeV, F z , = 11.5 ± 2.2 MeV The average 
masses are about 2 MeV above the thresholds of both B*B and B*B*. Recently, Belle Q also reported 
the observation of these two Z h states in T(55) —> [BB* + c.c.]-7r, and T(55) —> B*B*n. The discovery of 
the Zb states has inspired many interesting theoretical discussions. For example, it is suggested that these 
states can be molecular states of the BB* + c.c. or B*B* meson pairs [Ml- They are also proposed to be 
candidates of tetraquark states 0. In Ref. IE 15] the threshold enhancements are considered to be caused 
by cusp effects. 

Although the masses of these Zb states determined from the experimental fits are slightly above the thresh¬ 
olds, one should note that the masses are extracted by the Breit-Wigner parametrization. As emphasized 
in [9|. [l2j , if an S- wave shallow bound state exists below the threshold, the amplitude should not be param¬ 
eterized in the Breit-Wigner form. Using the line shape for a pure bound state, Ref. Q shows that the data 
on T(5S') —> hb(mP)n + n~ are consistent with the bound state nature of Z b K Furthermore, the observed 
enhancements in T(5S) —> [BB* + c.c]n and T(5S) —> B*B*n by Belle are very close to the thresholds 
of the £?(*)!?(*) systems. It is also found that the masses of the Zb states can be below the corresponding 
thresholds if these masses are extracted from data on T(55) —> B^B^tt Q. 

As a fact of observations, Zb states and their analogues in the charmonium sector ^,,(3900) [l6l - (l8l |. 
Z c (4020/4025) [HSJ and also the famous X(3872) appear to be strongly correlated to the thresholds of 
either B^ or £)(*)' pairs. This feature makes it natural to interpret these states as molecules. However, as 
was pointed out in Ref. [2l], [22j , it is difficult to understand the large production rates of these states in 
B-factories, e.g. X(3872), if these states are assumed to be loosely bound molecular states. In particular, 
the recent LHCb measurement of the ratio R^ = B ^x( 3 & 72 ytJf^ry) = ^ 0-64 ± 0.29 U seems not 

to support a pure D*°D° molecular interpretation of X(3872), since i?^ 7 is predicted to be rather small 
for a pure D*°D° molecule [2~dt . Meanwhile, a compact component inside such states can compromise both 
threshold phenomena and sizeable production rates. It is shown in Ref. j25|, [2(| that the radiative decays 
of X(3872) are not only sensitive to long-range parts but also to short-range parts of the wave function. 
The search for a hidden-beauty counterpart of X(3872), which is usually denoted as X,, is important for 
understanding the structure of X(3872). An effective field theory study shows that if X(3872) is a molecular 
bound state of D*° and D° mesons, the heavy-quark symmetry requires the existence of molecular bound 
state Xb of B*°B° with mass of 10604 MeV [27|. However, there is no significant signal of Xb near the 



threshold of B*°B° in Xb —> tt + tt~T(1S) [28] and in Xb —> wY(lS) j29|. Ref. [30] suggests that A'& may 
be close in mass to the bottomonium state Xbi{3P) and mixes with it. Therefore, the experiments which 
reported observing Xbi(3P) might have actually discovered Xb. 

Obviously, more experimental data and theoretical development are required to clarify the nature of these 
near threshold states. In Ref. 0 an effective field theo ry ( EFT) approach is proposed for the study of near 
threshold states (see also an independent study in Ref. [32]). In this framework the compositeness theorem 
can be incorporated with a determination of parameter Z which is the probability of finding an elementary 
component in the bound state, and the nature of near threshold states can be described by the presence of 
both molecular and compact components in their wave functions. 

The main purpose of this work is to study structure of Zb states by doing a combined analysis of data 
on T(5S) — > hb(mP)ir + ir~ and T(5S') — > B^B^n within EFT approach proposed in 0- Our work is 
organized as follows: in Sec. II, we recall the EFT approach proposed in Ref. [3l|. In Sec. Ill, we present the 
analysis of the T(55) —> hb(mP)n + TT~ transitions and in Sec. IV, the T(5S) —>• B^B^n. Our numerical 
results are presented in Sec. V. Finally, a brief summary is given in Sec. VI. 



FIG. 1: Feynman diagrams for the two particle scattering. The double lines denote the bare state. 


II. COMPOSITENESS THEOREM IN EFT 


In Ref. Jp|, we have developed an EFT approach which incorporates Weinberg’s compositeness theo¬ 
rem [33i. [341 . Here we recall some of the main points; more details can be found in Ref. [3lj| . Consider a 
bare state | B) with bare mass — Bq and coupling go to the two-particle state, where the bare mass is defined 
relative to the two-particle threshold. The two particles have masses mi, m 2 respectively. If | B) is near 
the two-particle threshold, then the leading two-particle scattering amplitude can be obtained by summing 
the Feynman diagrams in Fig. [1] Near threshold, the momenta of these two particles are non-relativistic. 
Therefore, the loop integral in Fig. [I] can be done in the same way as that in Ref. [35|, [Hj]. The loop integral 
can be written as 


Z = 


d D e 


(2tt) d )>o _ p/(2m{) +ie E — i° — P/(2m2) + ie 
f d D ~ x i _ i 

J {2tt) d ~ 1 e -l 2 /{2g) + ie 

-i2^T { 3 D )(Ait) L ^ (-2g,E - ie)^, 
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where /r is the reduced mass of the two particles, and E is the kinematic energy of the two-particle system. 
Obviously, the above integral does not diverge in D = 4. Using the minimal subtraction(MS) scheme which 
subtracts the 1/(D — 4) pole before taking the D —> 4 limit, one obtains 

Z = i^-{-2nE -ie) 1/2 . (2) 

Z7T 

It is interesting to note that, with the MS scheme no counter term is needed in the renormalization. We 
then have the two body elastic scattering amplitude for Fig. [T| 


A = - 
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If a bound state exists, we can have the following relations 
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where B is the binding energy, and Z is the probability of finding an elementary state in the physical bound 
state. Note that for the bound state, we mean a below threshold pole in the physical sheet. With Eq. ©, 
Eq. m can be re-expressed as 


where 


A= — 


E + B + £(£)' 


(5) 


t(E) = - s H± s /-2 fl E-u + - B)\. 

We can also express Eq. © in the form 

iA = ig 0 ■ G{E) ■ ig 0 , 

where G{E) is the complete propagator for the S -wave near threshold state 

iZ 


G(E) = 


E + B + t{E) +iT/2' 


( 6 ) 

(7) 

( 8 ) 


We have added a constant width T in the propagator, which can simulate the decay channels other than 
the bottom and anti-bottom mesons. From Eq. ©, one can find that the Feynman rule for the coupling 
between the near threshold state and its two-particle component is igo. Treating the binding momentum 
7 = (2 gB) 1 / 2 and the three-momentum of the two-particle state p as small scales, i.e., y,p ~ 0(p), one can 
then find that the leading amplitude Eq. ([5]) is at the order of 0(p~ 1 ). 
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FIG. 2: Feynman diagrams for Y(5 S) —¥ Z^iz —> hb(mP)mv, where the Zb states are produced in direct production 
processes. Solid lines in the loop represent bottom and anti-bottom mesons. 


III. T(5S) DECAYS TO h b (lP,2P)7r+7r- 

In this section, we study the decay T(5S) — > Z^tt — > hb(mP)mr in the EFT approach. Generally, in 

the decay T(5S I ) -A Z^ir, Zb states can be produced through both direct and indirect processes. In direct 
production processes, Zb states are produced directly via its compact component, while in indirect production 
processes a bottom and anti-bottom meson pair is produced first in the T(5S I ) decay and then rescatters 
to zjp. Similarly, the decay Zf ' 1 — > hb(mP)ir can proceed through both direct and indirect processes. In 
direct decay, Zb state will decay to hb(mP)TT directly. In indirect decay, zjp will first decay into a bottom 
and anti-bottom meson pair and then the meson pair rescatters into hb(mP)n. 
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The three-momenta of heavy mesons in decay T(5S') —> Z^n —> hb(mP)irir are small compared with 
their masses. Therefore, these heavy mesons can be treated as non-relativistic, and one can set up a power 
counting in terms of the small three-momentum p [3li . [dTj - liibt . From the power counting, one can find that 
if Zf' 1 contains a compact component, its production will be driven by this compact component [3l| (see 
also Ref. [H], E2|). In Figj2] we show Feynman diagrams where the Zb states are produced via compact 
components and decay through both direct and indirect processes. 

If zjp is a pure molecular state, its production should via indirect processes. The leading Feynman 

diagrams for indirect production of Z b ^ are shown in Fig. [3] Note that there are two kinds of indirect 
production mechanisms for Zb states. In Fig.[3fa,b), T(5S') first decays to a bottom and anti-bottom meson 
pair and pion in the same vertex, then the bottom and anti-bottom meson pair rescatters to Z b \ While 
in Fig. [3](c,d), T(5S') first decays to a bottom and anti-bottom meson pair, and after emitting one pion, 
the bottom and anti-bottom meson pair rescatters to zj ' 1 . It is shown in Ref. [l2| that both mechanisms 
contribute at leading order for the indirect production of Zb states. 
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FIG. 3: Feynman diagrams for T(5S) —> Z^n hb(mP)mv, where the Zb states are produced in indirect production 
processes. Solid lines in the loops represent bottom and anti-bottom mesons. 


As we are only interested in low energy physics, it is convenient to collect B mesons in a 2 x 2 matrix ['Kill 11 
H a = P:-d+P a , H a = -p* a -d + p a , P o W = (rW-, bW°), (9) 


where cr l are the Pauli matrices, and a is the light flavor index, P* and P a annihilate the vector and 
pseudoscalar heavy mesons respectively, and Pa annihilates the corresponding anti-particle. The leading 
effective Lagrangian describing the coupling of Z b states to the bottom and anti-bottom mesons can be 
written as that in Ref. Q 
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Z b HH = 
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2p2 


Tr [Z^H^Hb] + ^=Tr[(Z T Y ba HyHl], Z ab = 
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ab 


where Z ab annihilates Z ab , Z\ b creates Z ab , and go is defined in Eq. Q. The Lagrangian for the coupling of 
the P -wave quarkonia and the B mesons reads [371 ] 


C hb HH = y Tr [hl l H a a l H a ] + H.c. 


( 11 ) 
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The chiral Lagrangian for the B mesons and the S- wave quarkonia can be written as 0 

AfflxPT = ff T Tr[Hta i ff 6 ]4 6 - 5T Tt[fftH 6 a < ]^ + i 5 iTr[THtfft]A2 6 + ii^Tt[ T H>.^t ] + H.c. ) (12) 

where A d B = A(dB) — ( dA)B , T is the 2x2 matrix field defined as T = T(5S') • a + ih,(5S), and A^ is 
the axial vector pion current which is given by 

i , , iM ( -^7r° tt + \ 

A fl = ^d^-^) = -d fl M/F n + ---, £ = e*v , M=f^_ J , F« = 132 MeV. (13) 

We set g^ = 0.25 as in [§|, §2|. Note that our convention is different from that in Q, because a factor of 
\[2M has been absorbed into the field operator of the heavy meson in our convention (3lj| . then our g n is 
half of the value which is used in Ref. Q. The leading effective Lagrangian describing the Z^h^n interactions 
reads 


C Zbhb , = g z e^ k Zi b hl j A k ab + H.c., (14) 

which describes the direct decay of zj' 1 —> hb(mP)ir. 

Finally, we come to the vertex describing decay of T(5S') into Z b \. The corresponding Lagrangian to 
the leading order of the chiral expansion is given by [9j 


C rZb „ = grT(5S)ZtlA° ab + K.c.. 
Similar to Ref. [§j, we use the same coupling gr, g z for Zt and Z b . 


(15) 



FIG. 4: Feynman diagrams for non-resonant processes T(5 S) —> h b (mP)nn. Solid lines in the loop represent bottom 
and anti-bottom mesons. 


With the above effective Lagrangians and Eq. (J5J as the propagator of Z b \ one can then write out the 
amplitudes for all the Feynman diagrams in Fig.[2]and[3] We treat the loop integrals as was done in Ref. [37]. 
We present the relevant one loop three-point functions in Appendix A, and give all the amplitudes of Fig. [2] 
and [3] in Appendix B. In the following we address several points before ending this section. 

• As in Ref. 3, we assume that Zb only couples to BB* while Z b only couples to B*B*. We then 
find that there is a relative minus sign between iA43 a ,3b,3c,3d for T(5S) —> Z b n~ —> hb(mP)TT + TT~ and 
those for T(5S') — > Z' b + n~ —> hb(7nP)n + TT~ . It should not be surprising to find this relative minus sign, 
since if one assumes Zb(Z b ) couples to B*B*(BB*) with the same strength as that of Zb(Z b ) couples 
to BB*(B*B*), one would find that the meson loop amplitudes would be suppressed in heavy-quark 
spin symmetry world as noticed in [37]. 
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• Assuming that Z\, and Z' b are spin partners of each other, we can use the same Z for Zb and Z' b . In 
this way, we can reduce the number of free parameters in our fitting. 

• We show the Feynman diagrams for non-resonant contributions to Y(5 S) -A hb{mP)irn in Fig. [U 
Ref. [11| shows that the non-resonant diagrams do not satisfy the two-cut condition near Y(5S) region. 
Hence their contributions will not be enhanced by the kinematic singularity. We do not include their 
contributions in the present work, since they are suppressed by the heavy-quark spin symmetry. The 
experimental fits also find no significant non-resonant contributions [l|, Q . 



FIG. 5: Feynman diagrams for T(5 S) -A Z b \ r -A B^n. Solid lines in the loop and in the final state represent 
bottom and anti-bottom mesons. 


IV. Y(5S) DECAYS TO B w B (,) 7T 

In this section, we will study the decay Y(5 S) -A B^B^n in EFT. For the previous study one may refer to 
Ref. [h|, where the Zb states are assumed to be molecules. Instead of fitting data directly, Ref. pjJ constrains 
some parameters using data on T(5S I ) -A B^B^*\ and it then calculates the differential distribution for 
T(5S) -A £(*)£?(*) 7T as a function of invariant mass of the .£?(*) i?(*) pair. In this work, we give the amplitudes 
for T(55) -A 7t in EFT and constrain parameters by fitting the data directly. 

Similar to T(5S I ) — > Z b \ -A hb(mP)mr, Zb states can be produced through both direct and indirect 
processes. The leading order Feynman diagrams for these two different production mechanisms are presented 
in Fig. [5] The Feynman diagrams for the non-resonant contributions are shown in Fig. [G] We give all the 
amplitudes for Fig. [5] and [6] in Appendix C. 



(b) 



FIG. 6: Feynman diagrams of the non-resonant contribution to T(5S) -A B^B^x. The solid lines represent the 
bottom and anti-bottom mesons. 
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V. NUMERICAL RESULTS 


With the amplitudes given in Appendix B and C, we do a combined fit to data on T(55) —» 
hb(mP)n + n~ [|l[ Q and T(5S) —> B^B^n Q. Data on T(5S') —> B^B*n are nonvanishing below 
the £?(*)_£?(*) thresholds, hence we have to convolve the invariant mass spectra with detector resolution 
function. Data on T(5S') —> hb(mP)Tr + TT~ are collected per 10 MeV, so the invariant mass spectra should 
be convolved with detector resolution function and integrated over 10 MeV histogram bin. The detector 
resolution is parameterized by a Gaussian function with energy resolution parameter a = 5.2 MeV for 
T(5S) —> hb(mP)TT + tt~ [I| and a = 6 MeV for T(5S) —> B^B^n Q. To compare different scenarios for 
the structure of Zb states, we do the fit with two alternative schemes: 

a. We assume that Zb states are pure molecular states, then we have to set Z = 0 in the fit. In this way, 
only the diagrams in Fig. [3ja,c),[5jb,c) and [S] give nonvanishing amplitudes. 

b. We assume that Zb states contain substantial compact components, i.e., tetraquark component. It is 
shown in Ref. [22| that the production rate of a molecular state is proportional to its wave function square 
at the origin |'F(0)| 2 . Because the wave function of the molecular component in a loosely bound state 
spreads far out in space, |\I/( 0)| 2 is quite small, then the production rate of zjp through the molecular 
component will be suppressed. Therefore, we further assume that Z^' 1 is mainly produced through the 
compact component, and we set g\ = gi — 0 in the fitting. It is worth mentioning that Ref. pll . fiTTT 
demonstrate that the production of a near threshold state (by which we mean a mixture of the compact 
component and molecular component) is driven by the compact component. On the other hand, the 
hadronic decays of Z^ into hb(mP)TT will mainly go through the molecular component. This can be 
found from the power counting analysis. We treat the binding momentum 7 , the three momentum of the 
bottom meson ps and the four momentum of pion p w as small scales, i.e., they are all at the order of 0(p). 
Note that in the non-relativistic effective field theory, the propagator of the heavy meson is at the order 
of 0(p~ 2 ), and the measure of the one loop integration is at the order of 0(p 5 ). One can then find that 
Fig. (2a) is at the order of 0{p ~ 1 / 2 ), while Fig. [21(b) is at the order of O(p 0 ). Thus, as a leading order 
study, we set g z = 0 and neglect the contribution from Fig. [21(b). Up to now, we have shown that while 
the production of zjp is driven by the compact component, its hadronic decays mainly go through the 
molecular component. It is interesting to note that similar features are adopted for V(3872) in Ref. [2lj . 
By setting gi = g -2 = gz = 0, one can find that only the diagrams Fig. [ 21 (a) and El(a) give nonvanishing 
amplitudes, and the number of the relevant free parameters in this scheme is the same as that in scheme 
(a) (see Table. [J). 

We then compare our fitting schemes with that used in Ref. Although scheme(b) and Ref. joj use the 
same decay mechanism for T(5S') -A Zbir —>■ hbirn as shown in Fig.^a), there are some differences between 
them. The main difference is that Ref. Q sets Z = 0, while in scheme(b) we let Z to be a free parameter 
which satisfies 0 < Z < 1. As shown explicitly in Appendix B, the amplitude for Fig. Ha) is zero by setting 
Z = 0. Physically, by setting Z = 0, one assumes the Zb states as pure molecular states which do not 
contain compact components, hence they cannot be produced through the compact components. Therefore, 
if one uses Fig. Ha) to describe the decay mechanism of T(55) -A Zb'K —> hb tttt, one cannot set Z = 0 as 
in Ref. [9jJ. The consistent treatment is to let Z to be a free parameter which satisfies 0 < Z < 1 . On the 
other hand, if one assumes to be a pure molecular state, i.e., Z = 0, one should note that it can only be 
produced through indirect process. Therefore, for the pure molecular scenario, one should use Fig. Ha, c), 
i.e., scheme(a), instead of Fig.H a ) to describe the decay mechanism of T(55) -A Zbn -A hbTnr. 

Now we come to discuss the applicability of EFT. In the decay zY ] -A /i{,( 2 P) 7 r, the momentum of the 
pion is around 300 ~ 400 MeV in the energy region of our concern, hence the pion can be treated as soft, 
and one would expect the EFT expansion can converge fast enough. But in zjp —> hb(lP)ir, the momentum 
of the pion is relatively large and around 600 ~ 700 MeV. Based on naive dimensional analysis, Ref. 0 
warns that the EFT expansion may not be good enough for decay z'f' 1 -a hb(lP)TT due to the relatively 
large pion momentum. However, the results from the complete loop calculations can be more complex than 
the naive dimensional analysis. One may refer to Ref. [44j for an example. Generally, it is complex to study 
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the convergence of the effective field theory, and reliable conclusions can only be achieved once the complete 
higher loop contributions are available. Since such a kind of study is beyond the scope of the this work, we 
take a more pragmatic approach with two options in the fit. 

1. We use data sets of T(55) — > hb(lP 7 2P)n + n~, T(5 S) BB*n and T(5 S) — > B*B*tt in our fit. 

2. We use data sets of Y(55) —► hb(2P)n + n ~, T(5 S) — > BB*tt and T(5 S) — > B*B*tt in our fit. 

We choose an individual normalization factor for each final state in the fit. In this way, we need not to 
fix values of gr and gu- We present all the fitted parameters in Table. |U and we show the fitting results of 
ht(la), fit(2a) and fit(lb) in FigOEl Note that the width T in Table. Q]is not the total width, but the width 
defined in Eq. ©. 


TABLE I: Parameters for four fits. 


Fit 

52/ffi 

B 

B' 

rz b 

T z / 

Z 

X 2 /d.o.f. 

la 

0.049(15) 

0.11(12)eV 

27(58)keV 

2(l)keV 

1.9(1.9)MeV 

0 

110/58 

2a 

0.0017(69) 

12(21)keV 

0.14(7)MeV 

0.12(8)MeV 

0.59(27)MeV 

0 

72/45 

lb 

- 

0.19(22)eV 

1.6(1.8)eV 

5.5(1.8)MeV 

7.8(2.2)MeV 

0.42(12) 

81/58 

2b 

- 

0.38(65)eV 

0.51(86)eV 

6.1(2.8)MeV 

4.6(2.2)MeV 

0.42(18) 

69/45 


We give some brief discussions as regards our fitting results as follows 

• It is found in the experimental fits that the relative phase between Zb and Z' b in the hb(mP)Tnr channel 
is 180° [l, :2j- la fitting scheme (a), the relative minus sign between iM^, a , 3 c for Zb and Z' b can account 
for this relative phase. However, one can not find such a relative phase in amplitudes which are used 
in scheme (b). In our fitting, we find that scheme (b) gives a good fit only if such a relative phase is 
included. This may be attributed to gr (defined in Eq. (1151) 1 which has a relative minus sign between 
Zb and Z' b . We note that a very recent paper Ref. [45[ proposed an explanation for this relative minus 
sign. 

• From the fitting results of fit(lb) and fit(2b), one can find that the fitted parameters in fit(lb) and 
fit(2b) are close to each other. This indicates that the fitting results in scheme (b) are not sensitive to 
data on T(55 l ) —> hb(lP)n + TT~. Whether this means that the effective field theory can be successfully 
applied in T(5S') -» hb(lP)TT + n~ needs to be further investigated. Nevertheless, our numerical results 
show that such a possibility exists. It is also interesting to find that in scheme (b) the fitted binding 
energy and the width of Z b are close to those of Z' b . This seems to be consistent with the heavy-quark 
spin symmetry. 

• With all data sets, scheme (lb) gives much better fitting quality than scheme (la). Unfortunately, 
if data on T(55) —> /if,(lP)7r + 7r _ are dropped, the two schemes give almost equal fitting qualities. 
In this sense, it seems too early to claim conclusively that Z b states contain substantial compact 
components. However, a substantial compact component in Z b ^ can explain its large production rates 
in experiments. In contrast, a pure molecular state with the tiny binding energy as determined in 
scheme (a) is not likely to have large production rates in T(5S) decays. 

• The binding energies of the Zb states from the fit are generally very small. If we fix B = 0.1 MeV, which 
is the case for X(3872), and Z = 0.4 in fit(lb), we get a fitting quality y 2 = 90, which is still acceptable 
and better than fit(la). The other fitting parameters are B' = 0.23(14) MeV, T z b = 6.5(9) MeV and 
Tz' b = 5.6(9) MeV. This result also seems to be consistent with the heavy-quark symmetry. 

• One can also analysis data on T(5S) —> T(nS)TT + n~ in the EFT approach. However, different from 

hb(mP)Tr + TT~ and 7 r, the non-resonant contribution in T(?rS')7r + 7r is significant. It is im¬ 

possible to consider the interference with the non-resonant contribution correctly in one-dimensional 
analysis. To analysis data on T(5S I ) —> T(nS)ir + 7T~, one needs to fit the two-dimensional Dalitz 
distribution, which is beyond the scope of the present manuscript. 


8 





VI. SUMMARY 


We have done a combined analysis of data on Y(55) —> hb(lP, 2 P)tt + tt~, Y(5 S) —> BB*n and Y(55) —> 
B*B* 7r within EFT approach. With a combined analysis, we determine the resonance parameters of Zb 
states in two scenarios. In one scenario we assume that Zb states are pure molecular states, while in the 
other one we assume that Zb states contain compact components. It is found that by assuming that Zb states 
contain substantial compact components, one can have a better description of all data than by pure molecular 
assumption. By fitting the invariant mass spectra of T(5S') —> hb(lP, 2 P)tt + tt~ and T(5S I ) —> B^B^n, we 
determine that the probability of finding a compact component in Z^ is about 40% . It is also interesting 
to note that the probability of finding a compact component in Z^' 1 could be close to that in X(3872) which 
are around 26% ~ 44% gi, E3- 



FIG. 7: Comparison of the invariant mass spectra of hb(lP)ir and hb(2P)n in fit(la), fit(2a),fit(lb) and the experi¬ 
ment. The dotted line is the result of fit (la). The dashed line is the result of fit (2a). The solid line is the result of 
fit(lb). Data are from p|. 




FIG. 8: Comparison of the invariant mass spectra of BB* and B*B* in fit(la), fit(2a), fit(lb) and the experiment. 
The dotted line is the result of fit (la). The dashed line is the result of fit (2a). The solid line is the result of fit (lb). 
Data are from g| and have had background subtracted. 
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APPENDIX A: ONE LOOP THREE POINT FUNCTIONS 


The three-point loop functions we will encounter are 
d d t 


= H>/ 


M 12 M 23 


27T yfa 


^ ( 1 ° - + *)(<" +1-.2 + Sr , - + 6.2 - h , - ME + «] 

1 . c' — c . . 2 a + c — c' . 

tan (- ) + tan (- ) 

2\J a(c — ze) 2 \] a(c! — a — ie) 


(16) 


/ (1) (mi,m 2 ,TO 3 ,g) 
M12TO3 


(\/c — ie — \/c' — a ~ ie) + ^ I(mi,m 2 ,m 3 ,q), 

27rg 2 2 ^ 23 < 1 2 


(17) 


where = rriirrij /(to* + rrij) are the reduced masses, 6 32 = mi + m 2 — M, & 23 = m 2 + m 3 + < 7 ° — M, and 


V" 

/^(mi, m 2 , m 3 , g) is defined as 


a — | ; - I q~, c — 2/xi 2 6i 2 , d — 2/r 23 6 23 + t—^q 2 . 

m 3 J m 3 


(18) 


q l I (1 \mi,m2, m 3 , ( 7 ) = (-*) 


d d t 


{2*)* (< o _ JL + ie)(f 0 + bl2 + £. - ie)[/« + t 12 - fca - + ie] ' 


For more details, one can refer to Ref. 


APPENDIX B: AMPLITUDES FOR T(5S) ->■ h b (lP,2P)r 


The amplitudes for T(5S') —> Z^n hb(mP)n + n in Fig. [2] and [3] read 


*A4 2a = 


2 V 2 Z (j q t (j — (j } L / f - 


1 


£ ijk q l e* J (h b )e k ( T) 


F 2 E + B + t(E)+iT Zb /2 

x [I ( m B *, , ms*, <?) + /(wb, tob* , tob* , <?)] • 


(19) 


iM 2b = —i 


.grgzE-n 


£ ijk q i e* j (h b )e k (T). 


*A4 3a — — 


Fn E + B + t(E) + iT Zb /2 

2‘g 2 gig^ghE^ g / n ,.tp_ • u /2 _1 


( 20 ) 


^-(-2/rS-ie) 1 

7U P U + B + E(U)+ir Zb /2 


FI 

TV 

x £ ljk q l e* j (h b )e k (T)[I(-,m B *, m B , m B ,, 9 ) + I(m B , m B *, wb* , <?)]■ 


( 21 ) 
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iM$b = i 


V%Z ggig z E n 
27r FI 


/i(—2 pE — ie) 


^V2. 


£7 + B + E(B)+irz 4 /2 


£ vk q i e *i(h b )e k (r). 


( 22 ) 


bAd 3c = -16 


9 2 92ghgl r 
FI 


(^I (1 \mB*,mB*,mB,p) + p 2 <5 fcm + 


^/ (1) (m B ,TOB,TO S .,p) - I {1) {m B * ,p)^j p k p m 


E + B + Y,(E) + iT Zb /2 


x e ljk q l e* J (h b )e m (T) [J(m s *, m B , m B *, <?) + / (m B , m B *, m B *, <?)] • 


(23) 


*-M 3d = i4V2 J92 j£ J * (^I {1 \mB*,mB*,m B ,p) + ,mB,m B * ,p)^J p 2 S km + 


^ (1> (mB,rofl,m B .,p) - / (1) (mB*,mB,mfl*,p)j 
xe i ^ k q i e*^ {h b )e m {Y). 


'JZ 


E + B + Y,(E)+iT Zb /2 


(24) 


E is the energy defined relative to the £?£?* threshold. B is the binding energy, g is defined in Eq. a. 

is the energy of 7r , p is the three-momentum of the 7r , and 9 is the three-momentum of the 7r + . 
p = ^ le re duced mass. Note that the terms proportional to p k p m in A43 C and M 3 d will 

disappear in the heavy-quark limit, i.e., ms = ms*. This indicates that in the heavy-quark limit, the D 
wave decay of Y(5S) —> Z b tt is forbidden. We neglect the terms proportional to p k p m in the fit, since they 
will be suppressed by the heavy-quark spin symmetry. 

The amplitudes for T(5S) —> Z' b + tt~ —> h b {m.P)TT + n~ in Fig. [2] and [3] read 


7 .Ad 2a = 


2y/2Zg' grg^ghE^ 


1 


Fl E + B' + t(E)+iT Z ’j2 

x [I(m B * , m B *, m B *, q) + I (m B *, m B *, m B , 9 )] ■ 


B k r,i c *j (h b )e k (T) 


9 e 


(25) 


iM 


2b = -1- 




F 2 E + B' + t(E) + iT z J 2 


£ ijk q i 


:*3 


(h b )e k ( T). 


(26) 


iM 3a 


2g' 2 gig„g h E„ g _ 1/2 - J. - 

FI tt v E+B' + V(E)+iT z j2 

x £ ljk q l e* j (h b )e k (T)[I (m B *, m B *, m B *, q) + I (m B *, m B * ,m B ,q)\- 


(27) 


■ . . ■ ^J2Z g gig z E n . .] / 2 1 

*A4 3 f, = — z—-— g{—2pE — ie) ' -=- 

2tt FI m M ' £ + B' + E(£)+irW2 


£ v k q i e *i(h b )e k {T). (28) 


iA4 3c = 16 


9' 2 929hg 2 r 


E 2 


^/ (1) (mB*,mB*,mB*,p) + d (1) (m B ,m B * ,m B * ,p)) p a A fcm + 


(/“><■ 


mB*,mfl.,mB*,p) - / (1) (ms,mB*,mB*,|)J J p"p 


k rn 


1 


E + E' + E(E) + iT z //2 


x E ljk q' l e* : ’ (h b )e m (Y)[I ( tub * , tub * , ms * ,<?)+/ (mB* , ms *, tub , <?)] ■ 


(29) 
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iMzd = - * 4 \/ 2 ' 


9 929z9tt 
F l 


^/ (1) (?n_B*,TOB*,TOs*,p) + p*5 km + 

Vz 


(m B *, m B *, m B *, p) - j?V 

X£ iifc 9 i e *i( /ife ) e m(T). 


£ + B'+ e(e) + »rW2 


( 30 ) 


E is the energy defined relative to the B*B* threshold. B' is the binding energy, g' is the renormalized 
coupling constant which is defined in Eq. (0. We use g' here to distinguish from g , which is used in T(5S) —> 
Z b n~ -A hb(mP)TT + tt~ , since they may be different due to different binding energies. p is the reduced mass 
of the B*B* system. Other notations are the same with that in T(5S) — > Z b n~ -A hb(mP)n + n ~. We also 
neglect the terms proportional to p k p rn in the fit. 

In the above, we have assumed that Zb only couples to BB* while Z' b only couples to B*B*. We also 
assume that the probability of finding an elementary state in Zb and Z' b is the same. In other words, we use 
the same Z for Zb and Z' b . 


APPENDIX C: AMPLITUDES FOR T(5S) -a B (,) B W tt 


The amplitudes for T(5S) -A B + B*°tt in Fig. [5] and [6] read 


iM 5a 


9rgFn_ - VI - e ( T ) _ 

V2F„ E + B + X{E) + iT Zb /2 


(31) 


iM.5 b 


- <,V-yE - , e - (T) . e . (g . 0) 

y 2 ttF w E + B + T,(E) +iT Zb /2 


(32) 


iM5c = 4 g e* 3 (B*°)e m (T) (^I {1 \m B *, m B *,m B , <2w) + / (1) (to b *, m B ,m B *, qv)) qlS jm + 

^J (1) (m B ,mB,mB.,fe) - I (1) (m B .,mfl,m B .,9,)J gig” 


E + B + Y,(E) + iT Zb /2 


• (33) 


*M 6a = gi^e(T )-e*(B*°). 

r 


(34) 


iMfib = 


^4®! £ '(T) 6 fl(B*°) 

Ft 


E _ ^(^>5 - ^ ■ PbS ij ) + E 1 _ a {<iIp’b + 9tt • Pb? 3 ~ <&Pb) 


(35) 


E is the energy defined relative to the BB* threshold. E v is the pion energy, g is the three-momentum 
of the pion. p is the reduce mass of BB* system, pb and p B are the three-momentum of B + and B*°, 
respectively. A is the hyperfine splitting of the B mesons. 


The amplitudes for T(5S) -A B* + B*°tt in Fig. [5]and Fig. [6]read 


iMl5 a 


grg'En _ iyfz _ 

V2F v E + B' + t(E)+iT z j2 




(36) 
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iM 5b = - ff'V 2 P E ! £- £ '” e ‘(T )e* m (B*°)e* n (B*+). 

y 2ttF v E + B' + T,(E) + iT Z '/2 


( 37 ) 


iM.br. = *4 


g ,2 g2ff7r 

Ftt 


^/ (1) 11 (to b *, m B *, tob*, g7r) + 7 (1 ) (tob,tob*, m B *, g7r)) g^ fcm + 


»TlB*,mB*,mB*,g7r) _/(!)(. 
x^\ m (T)e* J (.B* 0 )e*’(£*+). 


tob , m B *, m B *, q*) I g£g™ 


E + B' + S(E) + ir z ,/2 


(38) 


iMlfta 


-ig 1 —£ imn e i {r)e* m {B* + )e* n (B*°). 

Ftt 


(39) 


LM 66 = -t4^e i (T)e* J '(B* + )e* fc (B*°) 

-T?r 


-Ett 


{-e^ m p^ql+S ikm p^qi) 


i ^ /A nmsrikrn^n ^kmnxij^m^n , „mik mi ~mki^mi , 'mil .:rri k ^ 

+ (E~ V 5 PbPb^ £ $ J P§Pb+ £ <ln<hr- £ QttFb+ £ Qtv PBj 


(40) 


E is the energy defined relative to the B*B* threshold. E n is the pion energy, is the three-momentum 
of pion. /i is the reduce mass of B*B* system, ps and pp are the three-momentum of B* + and B *°, 
respectively. A is the hyperfine splitting of the B mesons. 
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